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We solve analytically and numerically the generalized Einstein equations in scalar-tensor cos- 
mologies to obtain the evolution of dark energy and matter linear perturbations. We compare our 
results with the corresponding results for minimally coupled quintessence perturbations. Our results 
for natural (O(l)) values of parameters in the Lagrangian which lead to a background expansion 
similar to ACDM are summarized as follows: 1. Scalar- Tensor dark energy density perturbations 
are amplified by a factor of about 10 4 compared to minimally coupled quintessence perturbations 
on scales less than about 1000h~ Mpc (sub-Hubble scales). This amplification factor becomes even 
larger (> 10 6 ) for scales less than 100h _1 Mpc. On these scales dark energy perturbations constitute 
a fraction of about 10% compared to matter density perturbations. 2. Scalar- Tensor dark energy 
density perturbations are anti-correlated with matter linear perturbations on sub-Hubble scales. 
Thus clusters of galaxies are predicted to overlap with voids of dark energy. 3. This anti-correlation 
of matter with negative pressure perturbations induces a mild amplification of matter perturbations 
by about 10% on sub-Hubble scales. 4. The evolution of scalar field perturbations on sub-Hubble 
scales, is scale independent and therefore it corresponds to a vanishing effective speed of sound 
(c s jj = 0). It also involves large oscillations at early times induced by the amplified effective mass 
of the field. This mass amplification is due to the non-minimal coupling of the field to the Ricci 
curvature scalar and (therefore) to matter. No such oscillations are present in minimally coupled 
quintessence perturbations which are suppressed on sub-Hubble scales (c s $ = 1). We briefly discuss 
the observational implications of our results which may include predictions for galaxy and cluster 
halo profiles which are modified compared to ACDM . The observed properties of these profiles are 
known to be in some tension with the predictions of ACDM . 
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1. INTRODUCTION 

A wide range of cosmological observations indicate 
that the universe has entered a phase of accelerating ex- 
pansion. These observations include both direct geomet- 
ric probes of the expanding FRW metric and dynamical 
probes of the growth rate of matter perturbations. This 
growth depends on both the expansion rate and the grav- 
itational law on large scales. 

Geometric probes of the cosmic expansion include the 
following: 

• Type la supernovae (Snla) standard candles 0, 

• The angular location of the first peak in the CMB 
perturbations angular power spectrum^. This 
peak probes the integrated cosmic expansion rate 
using the last scattering horizon as a standard ruler. 

• Baryon acoustic oscillations of the matter density 
power spectrum. These oscillations also probe the 
integrated cosmic expansion rate on more recent 
redshifts using the last scattering horizon as a stan- 
dard ruler 

• Other less accurate standard candles (Gamma Ray 
Bursts @, HII star burst galaxies @)and standard 
rulers (cluster gas mass fraction [7j) as well as 
probes of the age of the universe [1] . 

Dynamical probes of the cosmic expansion and the grav- 
itational law on cosmological scales include: 



• X-Ray cluster growth data @ 

• Power spectrum of Ly-a forest at various redshift 
slices [lO, EH 

• Redshift distortion observed through the 
anisotropic pattern of galactic redshifts on 
cluster scales [ill, 113 

• Weak lensing surveys [T3l - il5j 

These cosmological observations converge on the fact 
that the simplest model describing well the cosmic ex- 
pansion rate is the one corresponding to a cosmological 
constant [l6j in a flat space namely: 

H{zf = Hi [O m(l + zf + (1 - O 0m )] (1-1) 

where H(z) is the Hubble expansion rate at redshift z, 
Hq = H{z = 0) and 51om the present matter density nor- 
malized to the present critical density for flatness. This 
form of H (z) and other similar, more complicated forms 
of it that are also consistent with cosmological data, are 
predicted by three broad classes of models: 

• Dark Energy Models attribute the observed ac- 
celerating expansion to an unknown energy com- 
ponent (dark energy [l~7l - [l9| ) with negative pressure 
and repulsive gravitational properties which dom- 
inates the universe at recent cosmological times. 
The simplest model of this class is based on dark 
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energy with constant energy density (in time and 
in space) the cosmological constant} la] 1 . This 
simple and physically motivated model however, 
is plagued by fine tuning problems [2l|. Alterna- 
tively, dark energy may be described by a mini- 
mally coupled scalar field (quintessence |2l| and k- 
essence[22|) which acquires negative pressure dur- 
ing an evolution dominated by potential energy 
Quintessence models also require the introduction 
of a fine tuned unnaturally small mass scale compa- 
rable to the Hubble scale today. The role of dark 
energy can also be played by the introduction of 
perfect fluids with proper equation of state param- 
eters [H| (eg Chaplygin gas [H, SI]) but this ap- 
proach is less motivated physically. 

• Modified Gravity Models are based on large 
scale modifications of General Relativity (GR) 
which lead to repulsive gravitational properties and 
accelerate the cosmic expansion. Representatives of 
this class of models include f(R) gravity 0, 
DGP models [28| and Scalar- Tensor (ST) theories 
29l l30j (for a historical review of ST theories see 
3l| ) . These models also require some degree of fine 
tuning of parameters but they are more appealing 
from a physical point of view. 

• Local Void Models assume the existence of an 
unusually large underdensity (void) [32l . |33| on 
scales of about 1 Gpc which induces an appar- 
ent isotropic accelerating expansion to observers 
located close to the center of it (within about 
20 Mpc). These models are in some sense minimal 
since they require no new theoretical input but they 
have two sources of fine tuning: the location of the 
observer at the center of the void and the statisti- 
cally improbable assumption of the formation of a 
1 Gpc void. 

Scalar- Tensor (ST) cosmological models H^] (extended 
quintessence [34j ) constitute a fairly generic representa- 
tive of modified gravity models. They are based on the 
promotion of Newton's constant to a non-minimally cou- 
pled to curvature scalar field whose dynamics is deter- 
mined by a potential U ($) and by the functional form of 
the non- minimal coupling F(<f>). 

The deviation of these models from GR is tightly con- 
strained locally by solar system observations and by 
small scale gravitational experiments [35|, H(|. These 
constraints however are significantly less stringent on 
cosmological scales j37l - l39l | and may also be evaded by 
chameleon type arguments [40l |. These arguments arc 
based on the fact that the effective mass of the non- 
minimally coupled scalar field depends on the local mat- 



ter density which is naturally high on Earth and in the 
solar system where most of the constraints are obtained. 
Such a high mass freezes the dynamics of the scalar field 
locally and allows it to mimic GR on solar system scales. 
On cosmological scales however, where the mean density 
is low, the field is allowed to dynamically evolve and drive 
the observed accelerated expansion. 

Another potential mechanism to relax the strong local 
constraints on ST theories involves the fact that GR is 
an attractor[4lj in the phase space of ST theories. An 
early cosmological deviation from GR is therefore allowed 
since it naturally leads to an evolution towards a late 
phase region that emulates GR in agreement with local 
space-time observations. 

The main advantages of ST cosmologies may be sum- 
marized as follows: 

• ST theories are natural alternatives to GR. Scalar 
partners to the graviton naturally arise in most at- 
tempts to quantize gravity or unify it with other 
interactions (eg the dilaton in supcrstring theories 
[42j or the radion in extra dimensional Kaluza- 
Klein[43[ and brane cosmologies (44) ) . 

• Even though ST theories respect most of the sym- 
metries of GR (local Lorentz invariance, conser- 
vation laws, weak equivalence principle) they are 
fairly general and several modified gravity theories 
may be viewed as effective special cases of ST theo- 
ries (eg f(R) gravity HH, Kaluza-Klein models (43j 
etc). 

• The direct coupling of scalar field to curvature 
(and thus to matter density) provides in principle 
a mechanism to evade the coincidence problem ie 
address the question 'Why did the accelerating ex- 
pansion start soon after matter domination in the 
universe f46j'. 

• ST theories in contrast to minimally coupled 
quintessence, natu rally allow superaccelerating ex- 
pansion rate [13, l48j . Superacceleration corre- 
sponds to an effective dark energy equation of state 
w c f[ < — 1. Thus, crossing of the phantom divide 
line w = — 1 is naturally allowed in ST theories. 
Superacceleration is equivalent to violation of the 
inequality 



dH(z) 
dz 



(1.2) 



1 Vacuum energy of quantum fields could in principle provide such 
a constant energy density[a| albeit with a much larger magni- 
tude than observed. 



In contrast, this inequality is enforced in most mod- 
els of GR with dark energy J4i| and corresponds to 
the null energy condition p c g > — p c g where /9 e ff, 
and Poff are the effective energy and pressure driv- 
ing the accelerated expansion [50| ■ 

• The growth of matter and scalar field perturbations 
is distinct in ST cosmologies compared to those 
based on GR. In particular, on sub-Hubble scales 
the growth of matter perturbations is driven by an 
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effective evolving gravitational constant [29|, |47| . 
The perturbations of the non-minimally coupled 
scalar field are scale independent and not negligible 
on small scales [2!| [5l| . This is in contrast with the 
case of minimally coupled quintessence where field 
perturbations scale as (5$ ~ H % ° an d are negli- 
gible on small sub- Hubble scales [HI, [53| (- ^> H). 
These amplified perturbations of non-minimally 
coupled ST dark energy (extended quintessence) 
trigger a local amplification in matter perturba- 
tions and could potentially provide a resolution to a 
puzzle of ACDM cosmology related to the concen- 
tration parameter of cluster and galaxy dark matter 
halo profiles [54[ . These profiles are observed to be 
significantly more concentrated than predicted by 
ACDM " 



r,r, 



1561 ] . Amplified perturbations of dark 
energy could naturally lead to a resolution of this 
puzzle [58^ . 

The growth of perturbations in ST theories and its quan- 
titative comparison with the corresponding growth in GR 
is the main focus of this study. This comparison can lead 
to the derivation of potential signatures of ST theories 
on the power spectrum and on other other observables 
related to the growth of density perturbations. 

Linear dark energy perturbations have been exten- 
sively studied especially in the context of GR [HI, HH 
and constrained by the CMB angular power spectrum 
spectrum (59l. [60j . In the context of minimally coupled 
quintessence it was found that the scalar field density 
perturbations exist on all scales but they are stro ngly 
scale dependent and negligible on sub-Hubble scales [53| ■ 
An interesting anti-correlation between dark matter and 
quintessence perturbations has also been pointed outj6l|. 
On larger (super-Hubble) scales, they can range up to 
about 10% compared to the matter density perturba- 
tions and they leave a trace on the low I multipoles of 
the CMB angular power spectrum through the ISW ef- 
fect [HI . The suppression of sub-Hubble scalar field den- 
sity fluctuations in minimally coupled quintessence origi- 
nates from the value of the non-adiabatic speed of sound 
in these models which is about equal to unity jfjjj, [HI . 
The sound speed determines the sound horizon of the 



dark energy fluid, l s = c s /H. On scales below this sound 
horizon, the perturbations vanish while on scales above 
l s dark energy can cluster. 

In ST theories, the non-minimal coupling of the scalar 
field to curvature perturbations (which in turn are driven 
by matter perturbations) leads to an amplification of the 
scalar field perturbations on sub-Hubble scales. Thus, it 
may be shown that on sub-Hubble scales the field per- 
turbations are scale independent [Hf and therefore 
the effective speed of sound c s <s> for ST field perturba- 
tions vanishes. As discussed below, the corresponding 
scalar field density perturbations are also amplified but 
they are anti- correlated with respect to matter perturba- 
tions ( (t't") < 0). In addition, as shown in the fol- 
lowing sections the ratio of the scalar field density per- 
turbations over the matter density perturbations on sub- 
Hubble scales is also independent of the scale and can 
become significant (up to about 10%) for cosmologically 
viable models. The goal of the present study is to demon- 
strate the above points using both qualitative analytical 
approximations and detailed numerical analysis. 

The structure of this paper is the following: In section 
II we derive the main equations that determine the evo- 
lution of the background and the linear perturbations in 
ST cosmologies. We also use analytical approximations 
to derive a few qualitative features of the solutions for 
the field and matter perturbations in these theories and 
compare them with the corresponding features of GR so- 
lutions. In section III we present a detailed numerical 
solution of the linear perturbation equations using a spe- 
cific form of the ST potentials which is able to produce 
an observationally viable expansion background similar 
to ACDM . Finally, in section IV we summarize our main 
results and discuss future prospects and extensions of this 
study. 



2. PERTURBATIONS IN SCALAR-TENSOR 
COSMOLOGIES 

We consider the following ST action in the physical 
Jordan frame |H, [H 



J 



S = 



R-Z&) g^d^dv® - 2C7($)) + S m [V 



(2.1) 



r 



where G is the bare gravitational constant, R is the 
scalar curvature of the metric g^ v and S m is the ac- 
tion of matter fields. In what follows we use units such 
that 8ttG = 1. The variation of the dimensionless func- 
tion F($) describes the variation of the effective gravi- 
tational constant. This variation (spatial or temporal) is 



severely constrained by solar system experiments [35Ll36| . 
The GR limit of ST theories is obtained either by fixing 
F(4>) = $ — 1 ($0 is a constant) or by freezing the 
dynamics of 3> using the function Z($) or the potential 
[/($). For example a large and steep Z($>) makes it very 
costly energetically for $ to develop a kinetic term while a 
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steep confining U(Q) (massive can make it very costly 
energetically for $ to develop potential energy. In both 
cases we have an effective freezing of the dynamics which 



reduces the ST theory to GR. 

Considering variation of the action (|2.1[) we obtain the 
dynamical equations 



J 



[ - ^uR 



+ Z($) ( a^a^ - -g^idop)' 



2Z(*)D* = -^i?-^(d Q $) 2 + 2^, 







r 



V M ^F($) - g^UFip) - g^U(^) , (2.2) 

(2.3) 
(2.4) 



where T^ v is the matter energy-momentum tensor T 1 *" = 
(2/v^g) x SSm/Sg^u. 

Considering a fiat cosmological FRW background 

ds 2 = -dt 2 + a 2 (t) [dr 2 + r 2 (d0 2 + sin 2 d<f> 2 )] (2.5) 

where matter is described by a prcssureless perfect fluid 
with density and pressure (p, p) = (p m , 0) we obtain the 
equations for the evolution of the background 

3F ■ H 2 = p m + ]f 2 -3HF + U = p tot (2.6) 

-2F-H = Pm + $ 2 +F - HF = Ptot + ptot (2.7) 

$ + 3H<I> = 3^(h + 2H 2 ) (2.8) 
d$ V J d$ 



p m + 3H p m = 



(2.9) 



where we have rescaled $ so that Z = 1 assuming that 
Z > 0. As discussed in detail in the next section it 
is straightforward to solve this system numerically with 
initial conditions corresponding to the time of recombina- 
tion (z ~ 1000) and find the evolution of the background 
homogeneous field $(f), scale factor a(t) and matter den- 
sity p m {t) ~ a(t)~ 3 . Notice that this is a system of three 
independent equations. In order to obtain them we sub- 
stitute H from Eq. (j2~7) in ([2T5|l and use Eqs. (j2~^l . 
(|2.8[) and (|2.9[) for the numerical analysis of the back- 
ground evolution. This substitution in (|2.8[) reveals a 



new effective mass for the field $ which depends on the 
background matter density p m (t). More details on the 
rcscaling and the initial conditions are provided in the 
next section. 

In order to obtain the evolution of perturbations we 
consider the perturbed FRW metric which in the Newto- 
nian gauge takes the form 

ds 2 = -(1 + 2<j))dt 2 + a 2 (l - 2i))5 ij dx i dx j (2.10) 



The linear gravitational potentials <f> and ip along with 
the scalar field perturbations 8$ are sufficient to fully 
determine the cosmological perturbations in these theo- 
ries. Using the perturbed metric (|2.10j) in the generalized 
Einstein Eqs. (|2.2[) it is straightforward to find the con- 
nection between the potentials 4> and ip 



4> = tp ■ 

§ • For F = 



(2.11) 

where F<& = For F = 1 (gR case) Eq. (I2~TT1) re- 
duces to the well known relation <f> = ip which is valid in 
GR in the absence of anisotropic stresses. It is straight- 
forward to obtain two additional differential equations 
for the perturbations (j>, ip and (5$ using Eq. (|2.10j) and 
perturbing 3> in (|2.2[) - (|2.3|) . These equations are of the 
form fil M, HI 



J 



£ + 2H£ + [ 3H - — ) - -(8p + 3Sq) = 



(2.12) 



+ 3H5<f> 



-5- + $* - -F*$ ) 5$ - - (2$ + 3H&U - - -^5R = 
a* ' 2 ' J 2 



(2.13) 



r 



where 



^3(H<f> + ^) 



(2.14) 



5 ^ = -2U + k ^ 



(2.15) 
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corresponds to the generalized Poisson equation, 

R = 6(H + 2H 2 ) (2.16) 
is the Ricci curvature scalar, 

fc 2 ' 



is the Ricci scalar perturbation and 



SR=2 



-£ - 4ff£ - 3H 



k 2 "i 



(2.17) 



J 



6q 



Sptot _ j_ 
F ~ F 

A2, 



$(5$ + \{F$R - 2£/*)<J$ + — + 2H^-{F^) + 

2 1 at 2 ai 

2 • 



2 fc 2 ' 

-3F 2 -ff+- — ) F,*(S$ 



- $ 2 </> - 2F0 - 2fl"F0 - -F£ 

3 



(2.18) 



is the effective pressure perturbation. Assuming specific 
forms for the field potentials F and U in Eqs. (|2.12[1 - 
(|2.13p . eliminating the gravitational potential <f> using 
Eq. (|2.1ip and using the background solution of the sys- 
tem (|2.6p - (|2.9|) . we may numerically solve for the pertur- 
bations ip(k,t), <5$(fc,i). After proper rescaling we solve 
the system in the next section with initial conditions at 
recombination (z ~ 1000) corresponding to an initially 
smooth scalar field $ and a properly rescaled gravita- 
tional potential ip in a background with small initial devi- 
ation from GR. Assuming that ip{k,t) and 5&(k,t) have 
been obtained by the numerical solution of the system 



(l2~T2>(t2~T3|) in the background solution of flU}, ([2TH]>. 
(|2.9j) . the gravitational potential 4>{k,t) can be obtained 
by using Eq. (|2.1ip that connects the two gravitational 
potentials in the context of ST theories. 

It is straightforward to obtain the density pertur- 
bations of both dark energy (<5$ = jf^) and matter 

(S m = in terms of the numerically derived perturba- 
tions <p(k,t), ip(k,t) and 5&(k,t) and the corresponding 
background. The total effective density perturbation 5p 
of Eq. (|2.15p may be expressed in terms of field and 
matter perturbations as [IT], [63| 



5fJL 



5 ptot Sp$ + 5pn 



1 

F 



F 



F 

F^R6$ 



U j9 6$-3H—(F, 9 5>f>) 



3H + 3iP 



+ (3HF - $ 2 )</>+ 



F 

+F£ + 5 Pm - -j^5$p r , 



(2.19) 



while the corresponding effective total density p is ex- 
pressed as (see also Eq. 



_ Plot P$+Pm 1 ( li 2 | TT Q TT Tp I 

"=y= F = f\2* + u - 3HF + pn 

(2.20) 

Using now the above discussed numerical solution 
for the perturbations <p(k,t), ip(k,t) and <5$(fc,i) and 
Eqs. (|2TT5]l . ([2~19| . (j2~2T)f we obtain the following den- 
sity perturbations 



from Eq. ([2~T9")l . 



Sp<s> = F ■ 6p,\s Pm =o 



Ptot = F ■ p 



(2.21) 
(2.22) 



fromEq. (l2~20l . 



5p$ _ Sp\s Pm=0 



/' 



(2.23) 



Ptot 

from Eqs. (|2~2"Tj) . (|2~22]) . 

<W = ^(tf/r P '-Sp:\sp m =o) = 5p to t~5ptot\s Pm =Q (2.24) 

where Sp^ is the total effective density perturbation 
as obtained from the Poisson Eq. (|2.15p and Sp\s Pm= o 
is obtained from Eq. ()2.19l) . The normalized matter 
density perturbation is obtained in terms of the field and 
metric perturbations as 

5 m = ^£HL = frg? - (2 . 2 5) 
Ptot M 



6 



where p is obtained from Eq. (|2.20p . 

Alternatively, the matter density perturbation, when 
normalized with respect to p m , i.e. S m = Sp m /p m , may 
be obtained on all scales by solving the differential equa- 
tion 



k f F 

-2H5 m + — U- -^<5$ ) -3(V>+2ffV) =0 ( 2 - 2 6) 



with appropriate initial conditions. 

Even though the system (|2.12j) . (|2.13p for the evolution 
of metric and field perturbations on a scale k can only 
be solved numerically, there are several useful qualitative 
conclusions that can be obtained by considering appro- 
priate limits of the scale k in these equations. There are 
four scales involved in Eqs. (|2.12j) , (|2.13p : 

• The physical scale - of the perturbations. 



where we have used <j> ~ ($ — <frj) H and where A is 
a proportionality factor necessary to fit the numer- 
ical solution. As expected, the above result implies 
that the scalar field perturbations are negligible in 
GR on sub-Hubble scales. 

The field density perturbations on these scales are 
obtained from Eq. (|2.19[) by setting Sp m = 
and considering only the scale dependent dominant 
term. We thus find 



tips 



k? 
«2 



k? 
a 2 



FF\ 



F 



F 2 



(2.30) 



The corresponding matter perturbation is obtained 
using Eqs. (|2~T5)) and (j2"30)) as 



Sp m = Sptot - <5p$ 



F 



F % 



(2.31) 



• The Hubble expansion rate H . 

1 /2 

• The mass scale of the potential U . 

• The shifted ST perturbation scale F 1 ^ 2 1 

In addition, in a cosmologically interesting setup, the ex- 
pansion of the universe is driving the time evolution of 
every physical quantity /. Therefore, |/| ~ H\f\. For 
scalar fields that can play a role in the present accelerat- 

1/2 

ing expansion of the universe we require U ^ ~ H . Thus 
each term in Eqs. (|2.12[) , (|2.13[) is determined by one of 
the three scales: |, H, F 1 ^ 2 ^. By identifying the terms 

that dominate in each range of perturbation scales - we 
may simplify the perturbation equations and obtain ap- 
proximate solutions for the corresponding range of scales. 
We consider the following ranges of perturbation scales: 



Sub-Hubble ST scales: 



F 



1/2 k 



> H, F* > 1. 



In this case the perturbation scale and the shifted 
ST perturbation scale are of the same order. Ignor- 
ing subdominant terms, Eq. (|2.13l) becomes 



The sub-Hubble ratio 



Sp.i 

Spr, 



is therefore obtained as 



Sp<s> 
Spm 



5$ 

Srr, 



F 2 



2F 



(2.32) 



This scale independence of the ratio j^- on sub- 
Hubble scales indicates that the effective speed of 
sound for ST dark energy perturbations is c s $ = 0. 
The fact that j^- < indicates an interesting au- 
tocorrelation between dark matter and dark energy 
perturbations which is also confirmed numerically 
in the next section. 

It is also straightforward to derive the equation for 
the evolution of matter density perturbations on 
sub-Hubble scales in ST theories. Eq. (|2.26|) on 
sub-Hubble scales takes the form 



2H5r, 



fc2 

«2 



F, 



-5$ 







(2.33) 



where we have (5$ given by Eq. (|2.28|) . Eliminating 
Kij) in Eq. ([233)) using Eq. (j23Tj) we obtain^ 



<5$ 



2V>)-F> 



(2.27) 



Using Eq. (|2.1ip to express </> in terms of ip we 
find 



6<f> 



FFc 



F + F% 



(2.28) 



Thus, on sub-Hubble scales the field perturbations 
5$ are independent of the scale k and can be a sig- 
nificant fraction of the total energy perturbations 
as demonstrated in the next section. The corre- 
sponding behavior in GR is very different. Setting 
F = 1 and ip = <f> in (|2.13D we obtain 



<y*~i4^-(*-*i)v 



(2.29) 



2H5 r , 



rS m 1 2F 



2 F2F 



3F% 



= 



(2.34) 



which provides the evolution of matter perturba- 
tions on sub-Hubble scales [47j . It may be shown 
that 



1 2F 



4F% 



F2F- 



where G ff is the effective gravitational constant 
in Cavendish-like experiments in the context of ST 
theories Therefore, Eq. ([234} has the antici- 
pated scale independent form (as in the case of GR) 
but Newton's constant G has been replaced by the 
effective ST gravitational constant G e s (remember 
that we have set 8irG — I). 
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• Sub-Hubblc GR scales: | > H > F^ 2 |, F<& < 1. 
If F$ <C 1, there is a range of sub- Hubble scales 
corresponding to 



H 2 a 2 
>* « -jfeT- « 1 



(2.36) 



for which the terms depending on the non-minimal 
coupling in Eq. (|2.13[) arc negligible compared 
to all other terms. For this range of sub-Hubble 
scales the scalar field perturbations are negligible, 
scale dependent and behave as in the case of GR 
(Eq. (J229])). On small enough scales however ie 



2„2 



H 2 a 



< F* <C 1 



(2.37) 



we reobtain the ST scale independent behavior of 
Eq. HQS) (see Fig. [HI lower pannel). 

• Super-Hubble scales: ^ <C H. In this case we may 
ignore the scale dependent terms in Eq. (I2.13[) 
to obtain S$ ~ ip. Clearly, there is no scale depen- 
dence for the perturbations on super-Hubble scales. 
Similarly, for matter perturbations we find (using 
Eq. (|2.26p ) S m ~ ip ie these perturbations are scale 
independent. 

The above qualitative features of the cosmological per- 
turbations will be confirmed and extended by the de- 
tailed numerical derivation of the perturbations that is 
presented in the next section. 



NUMERICAL SOLUTION 



these rescalings, the background and perturbation equa- 
tions are expressed in terms of the dimensionless barred 
quantities with no other changes (the form of all the equa- 
tions remains the same). In what follows we omit the bar 
(.7.) but we refer to the rescaled quantities. 

The next step is to define specific forms for the poten- 
tials F($) and U(&) to be used in the numerical solution. 
We consider the following form for F and U : 



ee (3+ (1 -jS)co£ 
U{$) = l + cxp(-A$) 



, Xf > 0.(3.8) 
(3.9) 



This form of F($) is consistent with solar system tests 
for $ ~ 0. Indeed solar system constraints of ST theories 
imply that [11] 



F 



< 10" 



(3.10) 



t=t 



where to refers to the present time. An additional advan- 
tage of this form of F($), J7($) is that they are similar to 
the corresponding functions that were reconstructed [48[ 
from the best fit expansion rate in the context of recent 
Snla data while F($) never becomes negative leading to 
instabilities. If $ remains close to during its evolution, 
as in the cases we consider, then we may expand around 
<I> = and keep terms up to $ 2 in Eq. (|3.8[) . Therefore, 
in our numerical analysis we use the simple generic form 



= 1 - A/* 2 



(3.11) 



The first step towards the numerical solution of the 
systems for the background evolution (Eqs. (|2.6I) - (|2.9[) .) 
and for the perturbations (Eqs. (|2.12[) . (|2.13[0 is a proper 
rcscaling to convert dimensional quantities to dimension- 
less ones. We use the following definitions of dimension- 
less (barred) quantities: 



t 


= H-}{ 




(3.1) 


H 


= H Ic fH 




(3.2) 


k 


= H re( k 




(3.3) 


v> 


= ipirf 




(3.4) 




= i(ii5$ 




(3.5) 


I'm 


~ 8nG Pm 




(3.6) 


u 


= H * ct u = 

8ttG 


= H 2 ci U 


(3.7) 



where we have set 8itG = 1 and H re f is a reference expan- 
sion rate used for the rescaling. Also ipi is the initial value 
(at recombination) of the metric perturbation ip. 2 Using 



2 The rescaling ip = ipiip allows to set ip = 1 in the initial condi- 
tions. Also note that ipi = ip(k, ti) is in general scale dependent. 



where A/ = (1 - 0)X f < 0. The forms of (([3~%]) and 

(|3.1ip ) and f/($) (|3 9[) for specific parameter values are 
shown are shown in Fig. [TJ Notice that for $ < 0.2 the 
two forms of F($) arc practically identical. 




FIG. 1: The potentials F($) and [/($) for the param- 

eter values A = 1, A/ = 5, p = 0.5 used in the numerical 
solution. In order to approximate F by Eq. (|3.1ip . in our 
numerical analysis we take the initial field value $i ~ 0.1. 
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We next solve the system ([21) ]t .(|!£g ]) .(|2T9" ll (after substi- 
tuting H in (|2.8p from (|2.6p ) to determine the back- 
ground evolution in the context of the above potentials. 
The initial time we consider corresponds to recombina- 
tion (a(ti) = j^j- = jfjQi) and is obtained by solving the 
Friedman equation d 2 = in the matter era as 



9ft 0m (l + ZiY 



1/2 



(3-12) 



where we have set i/ re f — Hq. The final time t/ of the 
solution is set to tf = 5 a value beyond the present to- 
The present time to is determined after the solution of the 
system by demanding that il m (to) = 3H(t°) 2 = ^ 0l] 



o.3 m) ^ m. 



Wc use the following initial conditions at tj correspond- 
ing to recombination 



4 



(3.13) 



obtained by assuming a matter era at early times 



*(*<) = 0.12 
*(ti) = 



a(U) = 



1 



1 + Zi 



1001 



(3.14) 
(3.15) 

(3.16) 



corresponding to initial conditions at recombination close 
to GR. We have checked that our results are robust with 
respect to reasonable changes of the above initial con- 
ditions. The parameters that need to be fixed for the 
solution of the background system (|2.6j) . (|2.8D . (|2.9[) are 
A/, fiom- In most solutions discussed in this section we 
set Oom = 0.3 and use the potentials (|3.9[) . (|3.11|) with 
A = 2, Xf = 5 (ST cosmology) or A/ =0 (minimally 
coupled quintessence). 

For the background system we find two solutions and 
select the one that is well behaved (the other diverges 
and is discarded). We also rescale the scale factor by its 
value at the present time (a(t) — > ^j) so that a (to) = 
1. This rescaling shifts somewhat the initial redshift Zi 
of our solution and the corresponding initial time given 
by Eq. (|3.12|l . We use the function t(z) obtained by 
solving numerically the equation = 1 + z, to find the 
value of the initial time at recombination corresponding 
to the rescaled solution a(i). This function is also used to 
convert plots vs time to plots vs redshift. Similarly, the 
inversion of the equation a = a(t) provides the function 
t(a) used for the conversion of plots of quantities vs time 
provided by the background solution to plots vs the scale 
factor a. 

We have verified that our solution reduces to ACDM 
for A = Xf = as expected. Deviations from ACDM 
are observed as we increase the value of A, thus steepen- 
ing the potential U and giving dynamics to Setting 



Xf = 0, we construct the numerical solution for the min- 
imally coupled quintessence field and the corresponding 
deviation from the ACDM expansion rate (Fig. [2]). For 
the potential in Eq. (|3.9[) we have checked that devia- 
tions from the ACDM expansion rate are less than 3.5% 
for any value of the parameter A with fixed value of ini- 
tial condition <&(t») — 0.1. At early times the scalar field 
dynamics is frozen due to cosmic friction, hence no devia- 
tion from ACDM arises. At late times, the constant part 
of the potential (|3.9[) prohibits a large deviation from the 
ACDM expansion rate despite of the somewhat increased 
kinetic energy of $. 



0.35 
0.30 
$ 0.25 
0.20 
0.15 



-3.0 -2.5 -2.0 -1.5 -1.0 
log,,, a 





FIG. 2: Background field dynamics for minimally coupled 
quintessence (top panel) and % difference between H and 
-Hacdm (bottom panel). Due to cosmic friction the field 
remains frozen at early times. Deviations from the ACDM 
expansion rate arise due to the field dynamics. During the 
evolution F$ = and F = 1. 



In Fig. [3] we plot the background dynamics for a non- 
minimally coupled quintessence field $ setting A/ = 5 
and A = 0.1. In ST gravity the scalar field is moving un- 
der the influence of the effective potential (see Eq. (|2.8|) , 
(EDS)) 

C/eff ($) = U{$) - ^ilF($) . (3.17) 

Assuming that F ~ 0(1) at all times this implies that at 
early times C/ C ff (*&) is dominated by the second term. The 
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effective mass of the field at early times is then m e g ~ 
(-RF^) 1 / 2 ~ ^/ 2 H. As a result, for \ f > 1 and ^ 
the scalar held is dynamically rapidly driven to its value 
corresponding to GR (<E> = 0) and remains performing 
oscillations of decreasing amplitude around the attractor 
$ = as the Universe expands. If $j = then the 
field remains at $ = (F ~ 1) until U e g($) ~ [/($). 
Therefore. <f> behaves as minimally coupled quintessence 
at early times and as non-minimally coupled quintessence 
at late times when the field is driven away from $ = 0. 
Since in the case shown in Fig. [3] we chose A = 0.1, the 
scalar 4> is not significantly driven away from 0, hence 
the present value of -F$ remains small enough to satisfy 

F 2 

the solar system constraints (Eq. (|3.10[1 ) since -j^-\t=t — 
3 x 10~ 4 . Thus, a small slope of the potential U($) is 
sufficient to secure that solar system constraints will be 
satisfied at present. 

Deviations with respect to the ACDM expansion rate 
can be significant at early times if the field begins far 
from <!> = 0. Therefore, large deviations from ACDM at 
early times can be avoided simply by tuning $^ so that 
~ 1. With such initial conditions, the system ap- 
proaches GR while the Hubble expansion rate becomes 
practically identical to ACDM at late times. When U 
starts dictating the field dynamics (late time evolution), 
such deviations are kept minimal by assuming a small 
slope of the potential £/($) (A = 0.1, Fig. [3]). As we in- 



crease the slope of the potential by increasing A (A = 2, 
Fig. 2]) we find small deviations (less than 1%) from the 
ACDM expansion rate at late times. Even for larger val- 
ues of A the deviations at late times are always below 3% 
but the dynamics of the scalar field $ is more interesting. 
In particular for steeper potentials U (A = 2) the dynam- 
ical evolution of $ leads to an amplified value of F\ at 
redshifts z < 0(1)- For example, in Fig. |Jwe plot the 
background dynamics with A = 2. In this case $ has sig- 
nificant late time evolution. This leads to an increased 



value of 



io- 



(Eq. 



violating solar system constraints 
3.10P but not cosmological constraints [38[ 



F % 



< 



0(1). 



(3.18) 



Thus, in the context of a chameleon mechanism the in- 
creased value of F\ has the potential of being consistent 
with observational constraints. Furthermore, and as dis- 
cussed in the previous section (Eq. (|2.32j) ). the value of 
F\ determines the significance of the dark energy density 
perturbations compared to those of matter. We therefore 
anticipate amplified dark energy perturbations (compared 
to GR) when the dynamics of $ is turned on by increas- 
ing the value of A. 
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FIG. 4: Background dynamics giving rise to the field contrasts depicted in Figs 6-11. From left to right and from top to bottom: 
1. Field oscillations about $ = 0; 2. % difference between H and Hacdm ; 3. Evolution of F($); 4. F% (~ 10 _1 at present). 
The parameters are A/ = 5 and A = 2. 



A remarkable feature distinguishing minimally coupled 
quintessence from the non-minimally coupled one is that 
the latter is able to cross the phantom divide line [3 [Hoi . 
l57l | corresponding to an effective dark energy equation 
of state w e s = — 1. In general, the effective equation of 
state for the scalar field <f> is given by (see Eqs. 
(EH)) 



I$2 



[/($)+ F + 2HF 
I$2 + „ 3H p 



(3.19) 



In Fig. [5] we show the evolution of the equation of state 
parameter for minimally (A/ = 0, solid line) and non- 
minimally (A/ = 2, dashed line) coupled quintessence, 
corresponding to GR and ST gravity respectively. At 
early times the background field oscillations give rise to 
divergences in the EOS parameter of $ but such diver- 
gences do not reflect on the Hubble expansion rate. At 
late times, the scalar potential £/(4>) becomes relevant 
for the field dynamics, the field starts growing and the 
EOS parameter oscillates around the phantom divide line 
W = —1. Crossing of the phantom divide line is allowed 
by all current cosmological observations and is in fact fa- 
vored by some of them (5(1 • This behavior is character- 
istic of ST gravities and cannot be achieved in minimally 
coupled quintessence [49) . 



-0.92 - 



-0.94 



-0.96 



-0.98 



-1.00 



-1.02 



General Relativity 



■ Phantom divide — 

s Scalar- Tensor gravity" 

1 2 

redshift z 



FIG. 5: Equation of state parameter (EOS) w$ for the back- 
ground field $ as obtained in General Relativity (A/ = 0, 
A = 2, solid line) and Scalar- Tensor gravity (A/ = 5, A — 2, 
dashed line). In Scalar- Tensor gravity the non-minimal cou- 
pling of the field results in the crossing of the phantom divide 
line (in red); for minimally coupled quintessence such crossing 
does not occur. 



We now use the background solution in Eqs. (|2.12l) . 
(|2.13[) to find the evolution of perturbations. We use ini- 
tial conditions corresponding to matter era in GR. Tak- 
ing also into account the rescalings p.4[) , (|3.5[) , the initial 
conditions used for the solution of (|2.12[) , (|2.13|) are of the 
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form 3 

6<S>(k,U) = 5<Z>(k,ti) = (3.20) 
= 1 ; V»(Mi) = (3.21) 

In order to find the value of the rescalcd wavcnumber k 
used in the perturbation Eqs. (|2.12p . (|2.13p correspond- 
ing to a given physical scale X p h _1 Mpc we express k as 
follows: 



2tt 



2tt 



3 x 10 5 kmsec 



X p h -1 Mpc H iei 



■H 



L Mpc lOOh kmsec -1 

(3.22) 

where c is the velocity of light and Hq = j^- is 
the present value of the Hubble parameter as pro- 
vided by the solution of the rescaled background system 

(&(&(& 

The evolution of the field (5<E> and metric perturbations 

ip vs the scale factor a (in logarithmic scale) for a per- 
turbation of wavelength A p = 30 h Mpc are shown in 
Figs. |B]and [3 In the top panel of Fig. [S]we plot the nu- 
merical solution for the evolution of 8$ as obtained in GR 
(solid line) and its predicted value on subhorizon scales 
(dashed line), given by Eq. (J223J). We set A ~ 25 in 
Eq. (|2.29[) in order to match the numerical solution. The 
bottom panel shows the numerical solution for the evo- 
lution of (5<E> in ST gravity (solid line) along with its pre- 
dicted value on subhorizon ST scales (dashed line), given 
by Eq. (|2.28p . As predicted by the analytical expression 
(|2~28|) . the field perturbat ion <5$ is a product of two oscil- 
lating modes: ip whose oscillations on sub-Hubble scales 
are driven by the term ^F^S® (see Eqs. ([2~T2")h (|27Tg|0 
with frequency 



1/2 



H 



and -F$ with frequency 



Xf-H 



(3.23) 



(3.24) 



This superposition of high and low frequency modes 
(aty 3> ^0 on sub-Hubble scales) manifests itself in Fig. 

m 
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FIG. 6: Evolution of the field perturbation 8& as obtained 
in General Relativity (top panel) and Scalar- Tensor gravity 
(bottom panel) for the scale X p = 30 h _1 Mpc. The numerical 
solution is the solid line and the analytical approximation on 
subhorizon ST and GR scales is the dashed lines. Oscillations 
in the background field $ induce oscillations in 5<E> through 
F,g> oc $. We use A/ = 5 for ST gravity and A = 2 in both 
cases. The spikes correspond to changes of sign of <5<F 




0.5 o.o 



FIG. 7: Evolution of the metric perturbation ip as obtained 
in General Relativity (A/ = 0, dashed line) and ST gravity 
(A/ = 5, solid line) for the scale \ v = 30 h _1 Mpc. We use 
A = 2 in both cases. 



3 Note that ip(k,ti) is actually ip(k,ti) = ^^'** j = 1 (rescaled 
by the initial gravitational potential). The initial gravitational 
potential corresponding to a scale invariant Harrison-Zeldovich 
spectrum obeys (j>(k,ti) 2 k 4 ~ S m (k,ti) 2 ~ k and therefore 
4>{k,ti) 2 ~ fc" 3 . 



In Fig. [7] we plot the evolution of the metric perturba- 
tion ip with the scale factor as obtained in GR (dashed 
line) and in ST gravity (solid line). Since ip is driven by 
terms having both frequencies uj^p and (see Eqs. (|2.12[) . 
(|2.18p ) we also see in the evolution of ip a superposition 
of high and low frequency modes. Notice that no such 
oscillations are present in the GR case (dashed line in 
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FIG. 8: Evolution of the density contrast 5<s> in ST gravity for 
a perturbation of wavelength X p — 30 h~ Mpc. The numer- 
ical solution is the solid line and the theoretical prediction, 
obtained using Eqs. (|2.28| l and (|2.30|l . is the dashed line. The 
sign of 5$ is negative, hence an underdensity in the DE fluid 
is created. We use the values Xf = 5 and A = 2. 



Fig. E). 

In Fig. [5] we plot the evolution of the density contrast 
<5<j> obtained in ST gravity with the scale factor. As ex- 
pected, the field contrast shows the same oscillations as 
(5$. Also expected is the good agreement between the 
numerical solution (solid line) and the analytical pre- 
diction (dashed line). The latter is obtained by using 
Eqs. (|2.23|) . (|2.30l) and the background solution for p to t 
(equations 12.201 12.22[) . On sub horizon ST scales we ob- 
tain a negative value of 5$, as opposed to the positive 
value of S m . This implies that the field and matter per- 
turbations are anti-correlated, i.e. an overdensity in the 
matter gives rise to an underdensity in the scalar field 
This result has also been obtained in the context of min- 
imally coupled quintessence [6l| . In that case however it 
is negligible on sub-Hubble scales. 



The evolution of matter overdensities is also affected 
by the introduction of a non-minimal coupling. This is 
demonstrated in Fig. [5] where we show the evolution of 
SmACDM in ACDM on a scale of 30h _1 Mpc (top panel) 
and the ratio x Sm — where S m corresponds to ST grav- 
ity (dashed line), minimally coupled quintessence in GR 
(solid line) and ACDM (dotted line). Clearly, the matter 
overdensity is amplified by about 10% in ST gravity while 
it is practically identical to ACDM in GR quintessence. 
This amplification of S m in ST gravity could be made 
even larger at the expense of introducing more signifi- 
cant variation of the background H{z) from ACDM than 
shown in Fig. [3l In Fig. [10] we show the correspond- 
ing scale dependence of the present matter overdensity 
S m ACDM{k 7 t ) as predicted by ACDM (top panel) and 
5m(Mo) — f or grp g rav jty (dashed line) and 
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GR quintessence (solid line). Clearly, the 10% amplifica- 
tion of matter perturbations in ST gravity is applicable 



FIG. 9: Top panel: Matter density perturbation 5 m ACDM 
in ACDM . Bottom panel: The ratio t — — where <5 m 

"mACLIM 

corresponds to ST gravity (dashed line), minimally coupled 
quintessence in GR (solid line) and ACDM (dotted line). 



on sub-Hubble scales while on larger scales the ampli- 
fication is negligible. This mild amplification of matter 
perturbations may be attributed to the corresponding 
amplification of dark energy perturbations in these the- 
ories which also affects matter perturbations despite the 
predicted anti-correlation. Indeed, the dark energy void 
in a cluster of galaxies reduces the negative pressure in- 
side the cluster and amplifies the gravitational collapse. 



In Fig. [TT] we plot the scale dependence of the ratio 
5<s>/Sm at present in GR (solid line) and in ST gravity 
(dashed line). Setting A/ = 5,A = 2, in the top panel 
we find that the GR solution (A/ =0, solid line) leads 
to negligible dark energy perturbations on sub-Hubble 
scales (<5$ ~ k~ 2 , Eq. (|2.29p ). In contrast, ST grav- 
ity (A/ = 2, dashed line) produces amplified, anticor- 
related with matter, dark energy perturbations on sub- 
Hubble scales. The ratio 4 s - is scale independent on 
these scales as predicted by Eq. (|2.32[) . In such mod- 
els where the dark energy perturbations can grow on all 
scales it may be shown [60( that the speed of sound c a 
vanishes. The anti-correlation is evident by the spike of 
the dashed line which corresponds to a change of sign 
of 5q> on sub-Hubble scales. For <C 1 (Fig. [3]) and 
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FIG. 10: Top panel: The scale dependence of the rescaled 
present matter overdensity 5 m ACDM(k,to)/ipi as predicted by 
ACDM . Note that tpi = ip(k, U) is in general scale dependent. 
Bottom panel: The ratio g to - j f° r ST gravity (dashed 

line) and GR quintessence (solid line). 

flat potential U (A = 0.1), there is a range of interme- 
diate sub-Hubble scales where ST perturbations behave 
as in the case of GR (5$ ~ fc -2 , lower panel of Fig. ITTj) . 
This case was discussed in section 3 using analytical argu- 
ments (sub-Hubble GR scales). On super- Hubble scales 
we obtain the anticipated scale independence. We have 
checked that the agreement between the numerical re- 
sult for the ratio (6$/6 m )o an d the analytical predictions 
on sub-Hubble ST (Eq. (|Z52jl ) and GR scales is very 
good. This is expected in view of the good agreement be- 
tween the numerical solution and analytical expressions 
5$ both in ST gravity, Eq. (g^ , and in GR, Eq. 



4. CONCLUSION 

We have investigated in detail, analytically and nu- 
merically, the evolution of dark energy and matter linear 
density perturbations in Scalar- Tensor (ST) cosmologies. 
We have found that the evolution of dark energy pertur- 
bations in ST cosmologies is significantly different from 
the corresponding evolution in minimally coupled (GR) 
quintessence. In particular, our results may be summa- 
rized as follows: 



FIG. 11: Scale dependence of the ratio 5<s>/8 m at present as 
obtained in GR (solid line) and ST gravity (dashed line). For 
the top panel we use A/ = 5 and A = 2, in which case there 
is no range of scales on which ST gravity behaves as GR. In 
order to display such an interval of scales more clearly, in the 
bottom panel we set A/ = 5 x 10 -4 and A = 0.1. 

• For natural (0(1)) values of the ST Lagrangian 
parameters which lead to a background expansion 
similar to ACDM , ST dark energy density per- 
turbations arc amplified by a factor of about 10 6 
compared to minimally coupled quintessence per- 
turbations on scales less than about 100h _1 Mpc 
(Fig. HU). 

• On sub-Hubble scales dark energy perturbations 
constitute a fixed fraction of about 10% compared 
to matter density perturbations (Fig. fTTj) . The 
fixed scale independent fraction implies that the 
effective speed of sound for ST dark energy is 
c s $ = 0. The corresponding fraction for minimally 
coupled quintessence perturbations scales as k~ 2 
and is about < 10 _4 % (Fig. ITTj) corresponding to 
c s $ = 1. 

• Scalar- Tensor dark energy density perturbations 
are anti-correlated with matter linear perturbations 
on sub-Hubble scales (Eq. and Fig. QT] 
where S m and 5$, are shown to have opposite signs). 
Thus clusters of galaxies overlap with voids of dark 
energy. 
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• The evolution of scalar field perturbations on sub- 
Hubble scales, is scale independent and involves 
large oscillations (Fig. [6]) induced by the ampli- 
fied effective mass of the field (Eq. (|5T7|0 . This 
mass amplification is due to the non-minimal cou- 
pling of the field to curvature and (therefore) to 
matter (Eqs. (|2.7|) . (|2.8|l ). No such oscillations are 
present in minimally coupled quintessence pertur- 
bations which are suppressed on sub-Hubble scales 
and vary as k' 2 (Eq. ([2~2"9"]) 1. 

• The evolution of matter density perturbations is 
affected by the introduction of non-minimal cou- 
pling (Figs. 9-10) and is amplified by about 10% 
in ST cosmology compared to minimally coupled 
quintessence and ACDM on sub-Hubble scales. 

• For small values of non-minimal coupling F t $ there 
is a range of sub-Hubble scales where the scalar 
field perturbations have a scale dependence similar 
to the case of GR (~ k~ 2 ). However, even in this 
case, for small enough scales ^ > -y- the field per- 
turbations become scale independent and enter the 
ST regime (Fig. [Til lower panel). 

These results have interesting observational conse- 
quences. In particular 

• Dark Matter Halo Profiles: ACDM predicts 
shallow low concentration density dark matter halo 
profiles for clusters and galaxies in contrast to ob- 
servations which indicate denser high concentration 
cluster haloes [1^]. The amplified anti-correlated 
with matter dark energy perturbation profiles can 
lead to a modification of the predicted by ACDM 
dark matter halo profiles. In particular, the dark 
energy voids in clusters of galaxies can amplify lo- 
cally dark matter clustering due to the local reduc- 
tion of negative pressure in the region of the cluster. 
A detailed investigation of this effect would require 
the solution of the full coupled nonlinear system for 
the evolution of dark energy and dark matter per- 
turbations under the assumption of spherical sym- 
metry. This is a straightforward generalization of 
the present study. 



• Large Scale Structure Power Spectrum 

P m {k): For a non-minimal coupling = 0(1) 

the ratio |f- ~ jwt) is scale independent for prac- 
tically all sub-Hubble scales (see Fig. [TT] upper 
panncl). Thus it would be hard to identify a scale 
dependent signature of dark energy perturbations 
on the matter power spectrum for such values of 
F^. For smaller values of the non-minimal cou- 
pling however, there will be a GR regime for large 
sub-Hubble scales where the dark energy perturba- 
tions are predicted to be scale dependent (Fig. [TT1 
lower panel) while on smaller scales we enter the 
ST regime where the ratio ^ becomes again 
scale independent. This transition from the GR 
regime on large sub-Hubble scales to the ST regime 
in small sub-Hubble scales may leave a trace (small 
glitch) on the matter power spectrum on a scale 

k ^ H 
a — pi/2 ■ 
■ * 

• Lensing by Galaxy Clusters: The lensing prop- 
erties of galaxy clusters may well be altered due 
to the presence of dark energy voids. It would be 
interesting to investigate the lensing signatures pre- 
dicted by the superposition of dark energy voids on 
galaxy clusters. 

In conclusion, the amplified and anti-correlated with 
matter, dark energy ST perturbations investigated in the 
present study provide a new direction of observational 
signatures for this class of modified gravity models. 

Numerical Analysis Files: The mathematica 
files used for the numerical analysis of this study 
and the production of the figures may be found at 
http://leandros.physics.uoi.gr/deperts/deperts.htm . 
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